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ABSTRACT 

We present evidence for very high gas fractions and extended molecular gas reservoirs in normal, near-infrared 
selected (BzK) galaxies at Z'-^l.S. Our results are based on multi-configuration CO[2-l] observations obtained 
at the IRAM Plateau de Bure Interferometer. All six star forming galaxies observed were detected at high 
significance. High spatial resolution observations resolve the CO emission in four of them, implying sizes 
of the gas reservoirs of order of 6-1 1 kpc and suggesting the presence of ordered rotation. The galaxies 
have UV morphologies consistent with clumpy, unstable disks, and UV sizes that are consistent with those 
measured in CO. The star formation efficiencies are homogeneously low within the sample and similar to those 
of local spirals - the resulting gas depletion times are ~ 0.5 Gyr, much higher than what is seen in high-z 
submm galaxies and quasars. The CO luminosities can be predicted to within 0.15 dex from the observed star 
formation rates and stellar masses, implying a tight correlation of the gas mass with these quantities. We use 
new dynamical models of clumpy disk galaxies to derive dynamical masses for our sample. These models are 
able to reproduce the peculiar spectral line shapes of the CO emission. After accounting for the stellar and dark 
matter masses we derive molecular gas reservoirs with masses of O.4-1.2x1O"M0. The implied conversion 
(CO luminosity-to-gas mass) factor is very high: aco = 3.6 ± 0.8, consistent with a Galactic conversion factor 
but four times higher than that of local ultra-luminous IR galaxies that is typically used for high-redshift objects. 
The gas mass in these galaxies is comparable to or larger than the stellar mass, and the gas accounts for an 
impressive 50-65% of the baryons within the galaxies' half light radii. We are thus witnessing truly gas- 
dominated galaxies at z ^ 1.5, a finding that explains the high specific SFRs observed for z > 1 galaxies. The 
BzK galaxies can be viewed as scaled-up versions of local disk galaxies, with low efficiency star formation 
taking place inside extended, low excitation gas disks. These galaxies are markedly different than local ULIRGs 
and high-z submm galaxies and quasars, where higher excitation and more compact gas is found. 
Subject headings: galaxies: formation — cosmology: observations — infrared: galaxies — galaxies: starbursts 
— galaxies: evolution 



1. INTRODUCTION 

Over the last decade, deep and wide multiwavelength 
galaxy surveys have been key in addressing critical issues 
regarding galaxy evolution. By using a variety of color se- 
lection techniques, and deriving photometric and/or spectro- 
scopic redshifts, different galaxy populations have now been 
probed to unprecedented detail up to at least z '--^ 3. Also, the 
history of cosmic star formation as well as the correspond- 
ing build-up of stellar mass has now been constrained using 
different observational approaches. 

Less well understood are the physical processes that regu- 
late the star formation activity in these young systems and the 
question of if and how galaxy star formation at high redshift 
differs from what is seen in the local Universe. While it has 
been established that star formation rates (SFRs) in galaxies 
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were on average higher in the past and that the cosmic SFR 
density peaked at z > 1 (see, e.g., LeBorgne et al. 2009), it 
is still a matter of debate how much of this increase is due 
to more frequent merging/interactions between galaxies, or 
other processes. Analysis of the ultraviolet (UV) and op- 
tical rest frame morphology and Ha velocity fields of high 
redshift galaxies suggest that the galaxies responsible for the 
bulk of the SFR density at z ~ 1 to 3 are disks (e.g.. Bell et 
al. 2005; Elbaz et al. 2007; Genzel et al. 2008). These find- 
ings confirmed early hints for the presence of extended disks 
at high redshift, based on damped Lya studies (e.g., Wolfe et 
al. 1986). Some authors, however, advocate a predominant 
role of mergers based on, e.g., irregularities in high signal to 
noise (S/N) 3D observations of emission lines (e.g., Flores et 
al. 2006). On the other hand, it is also found that distant 
star forming galaxies show irregular morphologies with high 
dumpiness, especially in the rest-frame UV (e.g., Cowie et al. 
1996; Chapman et al. 2003; Daddi et al. 2004; Elmegreen & 
Elmegreen 2005). In this context it is interesing to note that 
the existence of such bright clumps can explain the observed 
kinematic irregularities in high-z galaxies (i.e., without the 
need to require the presence of mergers, Immeli et al. 2004; 
Bournaud et al. 2008). More direct evidence for an impor- 
tant role of in situ, quiescent star formation in distant galaxies 
came from the estimate that the duty cycle of star formation 
is high in massive z '--^ 2 galaxies (Daddi et al. 2005; 2007a; 
see also Caputi et al. 2008), with typical durations of « 0.5- 
1 Gyr. This is much longer than the typical duration expected 
for starbursts ( <, 100 Myr) triggered by mergers, based on 
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simulations (e.g., Mihos & Hernquist 1996; Di Matteo et al. 
2008) and also on the observed gas consumption timescales 
for IR luminous merging galaxies (e.g., Downes & Solomon 
1998). 

A recent major step forward in characterizing the nature of 
star formation in distant galaxies was based on the discovery 
that star forming galaxies define a narrow locus in the stellar 
mass-SFR plane (in this paper we will refer to the ratio of SFR 
to stellar mass as the specific SFR, SSFR). This correlation 
has been estimated to have a slope of 0.7-1.0 in log space and 
a dispersion of a factor of only ^2 from redshifts z ~ 0.1 to 
z ^ 3.0 (Noeske et al. 2007; Elbaz et al. 2007; Daddi et al. 
2007a; Magdis et al. 2009). 

The normalization in this correlation, i.e. the SSFR at a 
given stellar mass, is higher by a factor of 30 at z = 2-3 com- 
pared to the local Universe. At even higher redshifts the situa- 
tion is unclear as the SFRs are typically less well constrained 
(as they are primarily based on UV luminosities). However, 
the available data is still consistent with the picture that the 
correlation extends to earlier epochs (3.5 < z < 7) with a nor- 
malization similar to that seen at z = 2 (Daddi et al. 2009; 
Stark et al. 2009, Gonzales et al. 2009). The existence of 
this correlation has major implications for our understanding 
of galaxy formation. For example, the cosmic evolution of 
the normalization of the relation (Daddi et al. 2009; Pannella 
et al. 2009; Gonzales et al. 2009) is similar to the cosmic 
evolution of the SFRD (Hopkins & Beacom 2006; 2008; Le 
Borgne et al. 2009, etc). It is possible in fact to describe most 
of the SFRD evolution as the combined result of the evolution 
of the stellar mass-SFR correlation and that of the space den- 
sities (e.g., the mass/luminosity functions) of galaxies. The 
possible peak at z '--^ 2.5-3 (Pannella et al. 2009; Magdis et 
al. 2009; Daddi et al. 2009a) of the SSFRs is thus likely more 
fundamentally related to the physics underlying galaxy for- 
mation than the peak at z > 1 of the SFRD evolution, as it is 
hinting at a change in the physical processes that regulate star 
formation. 

The remarkable tightness of this relation implies a very high 
degree of homogeneity in these galaxies. In other words, de- 
spite the physical complexity of the process by which stars are 
formed, for a normal, near-IR selected star forming galaxy, 
the star formation rate is ultimately linked to its stellar mass 
(which in itself is presumably related to the mass of the un- 
derlying dark matter halo). 

In order to make further progress in understanding how 
galaxy formation proceeds in distant galaxies, and to shed 
light on the origin of the stellar mass-SFR correlations, we 
need to start probing the properties of molecular gas reser- 
voirs in these galaxies, since it is the molecular gas out of 
which stars form. If the assembly of galaxies is mostly driven 
by mergers, one would expect most star formation in merger- 
driven starbursts with concentrated gas cores (like in submm 
galaxies) and very efficient star formation, and much less 
molecular gas and star formation in regular objects. If, on the 
other hand, most star formation is due to smooth gas flows, 
then star formation would predominantly take place in more 
quiescent disk galaxies, accreting gas from their surrounding 
medium. In addition, one could speculate that the tightness 
in the stellar mass-SFR correlation might be pointing to ho- 
mogeneous properties of the molecular gas reservoirs and star 
formation modes inside the galaxy populations. The increase 
in the normalization of this correlation with cosmic lookback 
time could then possibly be simply linked to an increase of 
the average gas fractions in galaxies, as well as in the ratio 



of the molecular (H2) over the atomic (HI) hydrogen fraction 
(see, e.g., Obreschkow & Rawlings 2009). The presence of 
massive gas reservoirs seem to be a major requirement for jus- 
tifying high SFRs, with the lack thereof possibly supporting 
instead the need for changes in the initial stellar mass function 
(IMF). 

Recently, we reported the detection of large molecular gas 
reservoirs in two near-IR selected massive galaxies at z ^ 1.5 
(Daddi et al. 2008). Based on the knowledge available at 
the time, i.e., the previously known empirical correlation be- 
tween CO line luminosity (Z-cq) ^'^'^ '■^'-^^ ^ luminosity (Lir) 
those galaxies should not have been detected in CO. How- 
ever, observations were still attempted based on the results of 
Daddi et al. (2005; 2007) that suggested a high duty cycle 
and long duration for star formation in massive high-z galax- 
ies. The subsequent detections of CO[2-l] emission indeed 
implied low star formation efficiencies {SFE, defined as the 
ratio Lih/Lqq). These SFEs are similar to what is seen in lo- 
cal spirals and much lower than previously seen in high red- 
shift galaxies. Subsequent studies revealed low-excitation CO 
[3-2] to [2-1] line ratios in these galaxies (Dannerbauer et al. 
2009) that are more typical of local spirals. 

Given that both galaxies were detected in CO, this early re- 
sult suggested commonplace large molecular gas reservoirs 
in distant massive galaxies (Daddi et al. 2008). However, im- 
proved statistics would clearly be needed in order to formally 
confirm the result and to estimate the variations in gas proper- 
ties among the high-z galaxy populations. Higher spatial res- 
olution observations are further required to constrain the sizes 
of the molecular emission; i.e., given the low SFEs one would 
expect extended reservoirs rather than the compact emission 
lines seen in submm selected galaxies (SMGs) for which Tac- 
coni et al. (2006; 2008) find typically half light radii of order 
2 kpc. 

In this paper we report on low-resolution (D-configuration) 
observations with the IRAM Plateau de Bure Interferome- 
ter (PdBI) of CO[2-l] (redshifted to 3mm) in four additional 
BzKs (near-IR selected galaxies at z 1 .5) thus tripling the 
number of sources in our sample. In addition, we report on 
higher resolution (~ 1-1.5"; B-configuration) observations of 
four galaxies. We also discuss the properties of these six 
galaxies with respect to the 'parent' BzK sample. We use the 
information gathered on these six galaxies in order to inde- 
pendently constrain the molecular gas mass and, conversely, 
the infamous CO luminosity to gas mass conversion factor 

"CO- 

The paper is organized as follows. In section|2]we present 
the new PdBI observations and their analysis. In Sect. |3]we 
discuss the relation of the six galaxies to the parent sample and 
the derivation of accurate SFRs and stellar masses. The SFEs 
are characterized in Sect.|4]and compared to other galaxy pop- 
ulations. We introduce numerical models of clumpy disks 
in Sect. |5] and use the models to calibrate dynamical masses 
for our sample in Sect. |6] This section also discusses impli- 
cations for the implied gas masses and the conversion fac- 
tors. A discussion and summary of our results is presented in 
Sect. I2] and |8] We apply a concordance WMAP3 cosmology 
throughout the paper. 

2. PDBI OBSERVATIONS AND RESULTS 

2.1. Target selection 

The six galaxies that were observed in CO were culled from 
the sample of Daddi et al. (2007ab) in the Great Observatories 
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Fig. 1 . — Velocity averaged maps of the C0[2-l] emission in our sample BzK galaxies. The (cleaned) images were obtained using natural weighting and are 
based on data from all observations (as indicated on top of the individual panels). Contour levels start at ±2cr and are in steps of Itr (see Tab.[T). The size and 
orientation of the beam is indicated in the bottom-left comer. The panels in the top two rows are 25" in size, those in the bottom row are 45" in size. The cross 
in each panel coiresponds to the VLA 1.4 GHz radio continuum position (see Tab.[T) and are ±2" (±17 kpc for z = 1.5) in size. 



Origins Deep Survey North field (GOODS-N). This sample is 
based on galaxies selected in the K-band down to A' < 20.5 
(Vega scale; ox K < 22.37 AB), from which we selected the 
population of star forming galaxies at 1.4 < z < 2.5 by using 
the "star forming" BzK color criterion of Daddi et al. (2004b) 
and requiring a detection in deep Spitzer imaging at 24/im. 



From these galaxies, we chose those that had a pre-determined 
spectroscopic redshift as required for follow-up observations 
in CO. Currently, 89 BzK galaxies of the Daddi et al. (2007) 
sample in GOODS-N have a spectroscopic redshift, from var- 
ious observing campaigns and from the literature. Among 
the six targeted galaxies, five redshifts were obtained through 
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Fig. 2. — ACS images in the F775W band (i-band) of the target BzK galaxies. C0[2-l] contours are the same as shown in Fig[T] The panels in the top two 
rows are 7.5" in size. The panels in the bottom are 15" in size. For reference, 1" subtends 8.5 kpc at z = 1.5. 



our GOODS-N campaigns at Keck using DEIMOS (Stern et 
al. in preparation). The redshift for BzK-12591 was de- 
rived by Cowie et al. (2004). This galaxy, having a strong 
bulge evident from HST imaging, has a possible detection of 
[NeV]A3426A emission line, suggesting the presence of an 
AGN. This galaxy is not detected in hard X-ray emission in 
Alexander et al. (2003) but exhibits some soft X-ray emission. 
In addition to the availability of a spectroscopic redshift, we 



also required the detection of the galaxy at 1 .4 GHz. This is 
important because the radio continuum allows us to derive an 
independent estimate of the SFR in the sources. 

2.2. Observations and data reductions 

Tab. [T] summarizes the PdBI observations of the six mas- 
sive, near-IR selected BzK galaxies at redshifts 1.4 < z < 1.6 
in the GOODS North field (D-configuration data for two 
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Fig. 3. — Spectra of CO[2-l] emission binned in steps of 25 km s"' . The yellow color indicates the regions where positive emission is detected. These regions 
have been used to derive total integrated fluxes and the velocity averaged maps in Fig.[T] The red lines show best fitting double Gaussian profiles to the spectra 
(see Tab. |2]for fit values). In the fits, we forced the FWHM of the two Gaussian functions to be the same within a given galaxy. For the four galaxies with 
the highest quality spectra we also show single Gaussian fits. Those yield significantly worse fits to the wings of the CO profiles, except for the face-on galaxy 
BzK- 16000 that is consistent with a single Gaussian. Velocities are computed with respect to the tuning frequencies of the observations as listed in Tab.[T] 



galaxies were already presented in Daddi et al. 2008). The 
frequencies during the observations were tuned to the ex- 
pected redshifted frequency of CO[2-l]. 

We used the compact D-configuration of the antennas for 
detection of the sources, providing the highest sensitivity and 



lowest spatial resolution (^5.5" at our typical frequency of 
90 GHz). Four of the sources were followed up at higher res- 
olution (^1.0-1 .5") using the more extended B-configuration. 
Typical integration times were 5-lOh per configuration. Most 
of the observations were obtained with six antennas, with the 
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Fig. 4. — Left: uv-coverage of the four target BzK galaxies that were observed also with the higher resolution B-configuration. Right: Signal amplitude 
versus baseline length in the uv plane. The solid curves show circular Gaussian models with FWHM as labeled. As a reminder, for a point/unresolved source the 
visibility amplitudes would be constant with uv distance. 



Table 1 
Observation summary table 



Source 


ra' 


DEC' 


Conf 


Obs. dates 


A inT 


Frequency 


Combined beam' 


rms 




J2000 


J2000 






h 


GHz 




/iJy* 


BzK-4171 


12:36:26.516 


62:08:35.35 


D 


04-08/2007 


8.2 


93.525 


2.05" X 1.55", PA = 64° 


41 








B 


01-03/2008 


11.2 








BzK-21000 


12:37:10.597 


62:22:34.60 


D 


04/2007 


6.3 


91.375 


1.87" X 1.60", PA = 64° 


51 








B 


01/2008 


6.2 








BzK- 16000 


12:36:30.120 


62:14:28.00 


D 


05/2008 


3.3 


91.375 


3.83"x3.22",PA = 47° 


54 








C 


1 1/2008 


2.8^ 














B 


01/2009 


4.4 














D 


05-07/2009 


6.5« 








BzK- 17999 


12:37:51.819 


62:15:20.16 


D 


05/2008 


6.0 


95.501 


1.52" X 1.30",PA = 96° 


41 








B 


01/2009 


6.4 








BzK-12591 


12:37:41.371 


62:12:51.06 


D 


05/2008 


5.4 


88.669 


6.37"x5.26",PA = 85° 


116 


BzK-25536 


12:37:28.357 


62:18:54.91 


D 


05/2008 


6.2 


93.754 


6.32"x5.10",PA = 90° 


72 



Notes: 

1: coordinates are from VLA 1.4 GHz continuum emission (Morrison et al. 2010). The VLA 1.4 GHz map has a resolution of 1.8", and the typical position 

accuracy for our sources is 0.20" 
2: equivalent six-antennas on source integration time. 
3: beam resulting after combining all available configuration and imaging with natural weighting. These are the beams sizes displayed in Fig.[T] 

4: noise per beam averaged over the full 1 GHz spectral range 
5: data taken with a primary beam attenuation of 24% while observing the source GNIO (Daddi et al. 2009b) 
6: data taken with a primary beam attenuation of 40% while observing another source 



exception of a few tracks in which only five antennas were 
available. The correlator has 8 independent units, each cov- 
ering 320 MHz (128 channels each with a width of 2.5 MHz) 
with a single polarization, resulting in a total bandwidth of 
about 1 GHz with both polarizations. 

We reduced the data with the GILDAS software packages 
CLIC and MAP, similarly to what is described in Daddi et al. 
(2008; 2009ab) and Dannerbaueretal. (2009). Bandpass cali- 
bration was performed using observations of the standard cal- 
ibrators J0418+380 and 3C273. Observations of J1044+719 
and J 11 50+497 were used for phase and amplitude calibra- 
tion. For flux calibration we used MWC 349 as a primary 
calibrator and 3C273 when the former was not available. A 



model was used to account for the frequency dependence and 
time variability of the emission of MWC 349. Given the typ- 
ical observed scatter between the observations and the model 
for MWC 349, and also comparing the results of flux calibra- 
tion to the implied antennas efficiency variations, we estimate 
that the flux calibration is accurate to 10% at 90 GHz. 

CO[2-l] emission has been successfully detected at the ex- 
pected frequency for all six galaxies. Fig. [T] shows the de- 
tection images (natural weighted, cleaned) averaged over the 
spectral channels where emission was detected. In all cases 
the peak flux is detected at S/N > 5. Fig. |2]show the same 
CO[2-l] contours overlayed on the Hubble Space Telescope 
(HST) images of the galaxies from the GOODS Advanced 
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Table 2 
Observed properties 



Source 


^Keck 


zco 


vco FWHM- 


, 1 

-"cop-ij 


S/N det' 


CO size 


CO blue/red'' 


PA blue/red S( 1.4 GHz) 








km s"' 


Jy km s"' 




kpc 


kpc 




MJy 


BzK-4171 


1.465 


1.46520 ±0.00058 


530 ± 32 


0.65 


8.3 


11.3 ±2.2 


1.7±1.7 




34 


BzK-21000 


1.523 


1.52133 ±0.00036 


444 ± 26 


0.64 


9.4 


6.2 ± 1.9 


5.3±1.3 


-25 ±13° 


43 


BzK-16000 


1.522 


1.52496 ±0.00028 


194 ± 18 


0.46 


8.7 


10.9 ±2.9 






19 


BzK-17999 


1.414 


1.41385 ±0.00032 


440 ± 33 


0.57 


9.4 


6.4 ± 1.9 


5.7±1.3 


-1±12° 


34 


BzK-12591 


1.600 


1.60019 ± 0.00076 


^60Q 


0.84 


5.3 








186 


BzK-25536 


1.459 


1.45911 ±0.00013 


^ 170 


0.36 


6.4 








28 



Notes: 

1: the accuracy of Keck redshift is generally better than 0.001, except for BzK-16000 and 21000 (see Section|23) 
2: derived from fitting double Gaussian profiles 



3: to determine the S/N of the detection we add in quadrature the S/N from point source fitting of the data in each configuration, using the spectral regions with 

positive line emission (yellow regions in Fig. [3) 
4: spatial separation between the red and blue component in the CO[2-l] spectra 



Table 3 
Rest-frame UV morphology 



Source 




n Sersic 


a/b 


PA 


FWHMcirc 




kpc 








kpc 


BzK-4171 


3.7 


1.0 


0.40 


86° 


3.4 


BzK-21000 


6.2 


0.5 


0.37 


-35° 


11.0 


BzK-16000 


5.1 


1.5 


0.97 




5.9 


BzK-17999 


4.7 


0.5 


0.42 


-39° 


7.7 


BzK- 12591 


4.5 


2.0 


0.85 


32° 


6.0 


BzK-25536 


3.0 


1.0 


0.32 


90° 


3.4 



Notes: 

The formal Galfit errors on the sizes and axis ratios coiTespond to an 
accuracy of better than 5%. The PA is measured East of North. FWHMcijt 
derived from the fit of a circular Gaussian profile. 




Camera for Surveys (ACS) campaign (Giavalisco et al. 2004). 
We note that in five of the six cases, the CO[2-l] is very 
well centered on the VLA 1 .4 GHz radio continuum positions 
(Morrison et al. 2010; crosses in Fig.|3]l and on the ACS rest 
frame ultraviolet images (Fig. |2]l. However, for object BzK- 
12591 the CO[2-l] detection is offset by about half a beam or 
^ 3". The CO[2-l] detection for this galaxy is the one with 
the lowest S/N ratio and is based only on low-resolution D- 
configuration data. We conclude that the offset is most likely 
due to a combination of the effects of noise (we expect uncer- 
tainties of 0.6-1" at Ict in the position) and phase instabilities 
(the interferometric seeing was about 1.9") during the late- 
spring observations (Tab.[T]i. 

2.3. C0[2-l ] detection images and spectra 

Fig. [3] shows the resulting CO[2-l] spectra binned with a 
resolution of 25 kms^^. The spectra were obtained by fitting 
the data in uv space with the following procedure: we first 
averaged the uv data over the spectral range in which CO has 
been detected (see the yellow shaded regions in Fig.O to de- 
termine the spatial centroid of the CO emission. Keeping this 
centroid fixed for all channels, we then extracted the spectra 
with 25 km s^^ binning by fitting each channel in the uv data. 
For those objects for which we also have B-configuration data 
and that appear resolved (BzK-4171, 21000, 16000 and 1799, 
see next section) we used a circular Gaussian model with a 
fixed FWHM (derived from the velocity averaged data). For 
the remaining two targets (BzK-12591 and 25536) we fitted 
point source models. 

To determine the properties of the CO[2-l] emission (ob- 



FlG. 5. — 3-color RGB imaging of the six CO detected BzK galaxies, based 
on ACS imaging of Giavalisco et al. (2004), release 2.0. Red is F850LP, 
green is F775W, blue is F435W. We have used a 2-pixel tophat smoothing to 
improve the source visibility. Images are 3" in size. 



served flux, velocity centroid and FWHM) we fitted Gaussian 
functions to the observed spectra (Fig.O. Several of the spec- 
tra show clear evidence for double-peaked emission profiles, 
suggestive of rotation. We obtained fits for all galaxies with a 
single Gaussian and with two Gaussian functions (in this case 
fixing the FWHM in each component to the same value). In 
most cases we find that a fit with two Gaussian profiles pro- 
vides a better description of the data. This includes objects 
like BzK-21000 and 17999 where, despite the fact that the ev- 
idence for a double peak profile is only tentative, a fit with 
a single Gaussian is disfavored as it is not able to accurately 
reproduce the rapid decrease of the flux at the edges of the 
spectra. In no case we do find evidence for continuum emis- 
sion. We have combined all line-free channels in all galaxies 
to derive an average continuum at 3.3mm of < 100/iJy at the 
3fT level. This is consistent with the expectation, given the 
moderate SFRs in these galaxies (see Section[3]i, that the con- 
tinuum emission should be at the level of 20-30/iJy. 

Tab. |2] summarizes the results of the spectral fitting. Typ- 
ical integrated CO[2-l] fluxes are 0.4-0.8 Jy km s^'. The 
fufl width at half maximum (FWHM) of the CO [2-1] emis- 
sion is typically in the range of 400-600 km s^', with two 
notable exceptions: BzK-16000 and 25536. For BzK-16000 
the narrower velocity range is most likely due to the fact that 
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BzK-4171: C0[2-l] pver r775W 
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BzK-1799S: CO(E-lJ over FYTSW 
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Fig. 6. — Contours of CO[2-l] are overlayed, separately for the blue (approaching) and red (receding) velocity ranges in the spectra, on ACS F775W images 
(5" size). Contours start at ±2cr and increase in steps of Icr. The displacement of the red and blue parts of CO[2-l] for BzK-21000 and 17999 appeal' to be well 
aligned with the rest-frame UV major axis (dashed lines). 



the galaxy is observed nearly face-on. For BzK-25536 this 
is likely, at least in part, due to the relatively low S/N ratio 
of this source (see later in Sect. |6l). The velocity centroid of 
the CO[2-l] emission can be used to derive a CO[2-l] red- 
shift for the six galaxies. The CO redshifts agree well with 
the Keck redshifts to within Az = (zco — ZKeck) = 0.001, the 
typical accuracies of the latter. Two objects show a larger off- 
set: BzK-16000 has Az = 0.003 that can be, at least partly, 
explained due to the lower accuracy of the Keck redshift. In 
this case, the optical redshift was determined through weak 
absorption lines in the UV rest frame (all other redshifts are 
from the detection of [OII]A3727 emission lines); BzK-21000 
has Az = —0.0017, but its [Oil] emission line is close to a 
very bright, telluric OH emission line. 

2.4. Estimation of the spatial size of C0[2-l ] 

When comparing the D-configuration and B -configuration 
data (their respective uv coverages are shown in Fig. |4]i, im- 
aged separately, we find brighter peak flux densities in the 
D-configuration only data. This already indicates that the 
sources are spatially resolved. In order to quantitatively con- 
firm this finding and to estimate the intrinsic size of the CO[2- 
1] emission we have fitted models to the uv data. We have 
averaged the data over the spectral ranges showing positive 
emission (yellow-shaded regions in Fig. [3]). Given the mod- 
est S/N ratio of the observations, we fitted circular Gaussian 
models to the data to limit the number of free parameters (four 
in this case: the spatial center RA-DEC positions, the total 
flux and the Gaussian FWHM). Fig. |4] (right) shows the total 
CO[2-l] amplitude binned in ranges of uv distances, covering 
baselines from ^ 50m to ^ 300-400m. The overplotted model 
corresponds to the expected behavior in uv space of circular 
Gaussian profiles with the fitted sizes. In all four cases we 
find that the sources are spatially resolved with a significance 
of > 3(7 (5f7 for BzK-4171). The FWHMs of the circular 
Gaussian functions range within 0.7-1.3", corresponding to 
6-11 kpc at the redshifts of the galaxies. 

2.5. Comparison to ACS imaging 

It is interesting to compare the CO[2-l] sizes to those from 
the HST imaging. In order to do that, we used the Galfit code 
(Peng et al. 2002) to characterize the ACS z-band images 
(F850LP filter), corresponding roughly to U-band rest frame. 
Tab.[3]summarizes the results of the Sersic profile fitting. Pre- 
liminary results for the Sersic profile fitting of BzK-21000 and 



4171 were presented also in Daddi et al. (2008), and are con- 
sistent with the new measurements based on HST images that 
have twice the integration time (ACS GOODS release 2.0). 
Most of the sources are consistent with nsersic ^1-5 profiles 
that are typical for spiral disks. For BzK-12591 and 16000 
we find higher Sersic indexes of 1.5 and 2. Inspection of the 
HST images (Fig. |5]l shows that these galaxies appear to be 
disk galaxies with a significant bulge component. We cannot 
exclude that bulges might exist in some of the other galaxies 
as well, and deep high resolution near-IR observations (e.g., 
with WFC3 on HST) would be necessary to check this. Two 
galaxies have a Sersic index lower than 1: BzK- 17999 and 
21000. For these two galaxies most of the UV light is con- 
centrated in bright clumps, implying a light distribution flat- 
ter than an exponential. This is consistent with the appearance 
expected for gas rich disks (Bournaud et al. 2008). We notice 
that in most cases the target BzK galaxies display a significant 
inclination in the ACS imaging, with axis ratios often below 
0.5, i.e. they are not circularly symmetric. This implies that 
the size estimate derived from Sersic model fitting might not 
be directly comparable with the cruder CO[2-l] sizes derived 
from (symmetric) circular Gaussian fitting. While it is not 
meaningful to fit our CO[2-l] observations with more com- 
plex spatial shapes due to the limited S/N ratios, we decided 
to fit also the ACS images with circular Gaussian models to 
test for consistency between CO[2-l] and rest frame UV sizes. 

The results of this analysis are listed in Tab. [3] For three 
galaxies the size measurements agree within l-2a. The ex- 
ception is BzK-4171 for which the CO[2-l] size estimate is 
about three times larger than the UV rest frame size, a differ- 
ence significant at the 3.4cr level. This is also the object with 
the most accurate CO[2-l] size determination in our sample. 
If we average over the sample in order to increase the S/N ra- 
tio, we find that the ratio of CO[2-l] to rest-frame UV sizes is 
1.24 for all the four galaxies and 0.96 when excluding BzK- 
4171. Therefore, it appears that CO[2-l] and rest-frame UV 
sizes are quite similar, within 20% on average. 

2.6. Evidence for rotation 

We have analyzed the high-resolution CO[2- 1 ] data of three 
of our galaxies (as one, BzK- 16000, appears to be seen face- 
on) to check for signatures of rotation. Solid evidence for 
galactic rotation in near-IR selected massive BzK galaxies has 
recently been derived using Ha observations from the SINS 
survey (e.g., Forster-Schreiber et al. 2009). Clearly, it would 
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Fig. 7. — Left panels: CO[2-l] images with superimposed pseudoslits that 
are used to extract position-velocity diagrams. Right panels: the respective 
position velocity diagrams. The arrow shows the spatial and velocity scales 
in the images (1" corresponds to 8.5 kpc at z = 1.5) . The data were binned 
in velocity (width: 100 km s"'). 



be important to confirm such findings based on the dynamics 
of the molecular gas that should trace the underlying velocity 
field. 

In Fig.|6]we show CO[2-l] contours imaged separately for 
the blue- and red-wing velocity ranges of BzK-4171, 21000 
and 17999. For BzK-21000 and 17999 we observe a signif- 



icant spatial displacement of the red (receding) versus blue 
(approaching) component. Their displacement is, not unex- 
pectedly, aligned with the galaxy major axis as determined 
from the UV rest frame observations. Also, the separation of 
the blue and red components is found to be quite consistent 
with the overall UV sizes of the galaxies. We conclude that 
we find evidence for galactic rotation in the disks in at least a 
couple of our sources. 

For a more direct assessment we obtained position-velocity 
diagrams for our galaxies (Fig.|7]i. In order to obtain sufficient 
S/N we have binned the data over 100 km s^' and extracted 
position-velocity diagrams using pseudo-slits at the orienta- 
tion defined by the major axis in the ACS imaging. As can be 
seen from Fig.|7] we again find evidence for rotation. In par- 
ticular, BzK-21000 shows the expected S-shaped structure in- 
dicative of disk rotation, consistent with a flat rotation curve. 
A less pronounced S-shaped rotation structure is also seen in 
BzK-17999. 

In the case of BzK-4 171 the approaching and receding parts 
of the emission are not significantly offset (Fig. |6]l, and no 
clear evidence for rotation is present in the position-velocity 
diagram. We explored position velocity extractions at differ- 
ent orientations but could not find evidence for rotation at any 
position angle. We recall that for BzK-4171 we also got the 
puzzling result that the CO[2-l] size was found to be substan- 
tially larger than the rest-frame UV size. This object has the 
most pronounced double horn profile in our sample, with a 
very significant central minimum reaching about 15% of the 
peak flux of the two symmetric horns. This galaxy is the most 
compact object in our sample in the rest-frame UV, and the 
half light radius from the Sersic fit is consistent within the 
limits with the spatial separation between the red and blue 
component. Higher resolution and higher S/N observations 
are required to constrain the rotation in the disk of BzK-4 171. 

3. GENERAL PROPERTIES OF THE SAMPLE 

In this section we discuss the overall properties of the CO- 
detected galaxies. We compare their colors, stellar mass and 
star formation rates to the full sample of near-IR selected 
galaxies at similar redshifts to investigate whether they are 
representative galaxies for their mass. We further discuss the 
accurate derivation of the star-formation rates through multi- 
wavelength indicators. We also explain the measurements of 
the stellar mass and stellar population parameters through fit- 
ting of synthetic models to the observed multicolor spectral 
energy distribution (SED). 

3.1. Comparison to the parent sample 

We have checked whether our source selection criteria 
(Sect. 12.11 ) imposed substantial biases on the choice of CO 
targets. Fig. [8] shows the stellar mass-SFR diagram for the 
CO-detected galaxies in comparison with the full sample of 
BzK galaxies in GOODS-South. In the following we use 
the GOODS-South sample for comparison even if the CO- 
detected galaxies are in the GOODS-N field as our GOODS- 
South catalog reaches 1.5 mag deeper in K, down io K <22 
Vega. In order to use the same estimators for all galaxies in 
this plot, we use the empirical relations by Daddi et al. (2004) 
to estimate stellar masses and SFRs from UV luminosities 
in all objects (we discuss more accurate derivations of these 
quantities for the CO-detected galaxies below). From this 
plot it appears that the targeted galaxies are massive, but that 
they lie right on top of the stellar mass-SFR correlation. The 
reason why they are massive systems is due to the fact that we 
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Fig. 8. — The location of the six BzK galaxies (large filled circles) in the 
stellar mass-SFR diagram. SFRs are derived from the UV luminosities at 
1500A, connected for dust reddening. We show typical eiTor bars in the SFRs 
of the CO detected galaxies based on the results in Sect|3] The smaller points 
(squares; filled symbols are for spectroscopic redshifts, which are 30% of 
the sample) show individual BzK galaxies from the GOODS South sample 
of Daddi et al. (2007ab) to probe smaller stellar masses (only star forming 
BzK galaxies, or sBzKs, are shown). The magenta line shows the derived 
correlation between stellar mass and SFR. Filled magenta squares show the 
SFRs measured from radio stacking as a function of stellar mass (see also 
Pannella et al. 2009). The light blue lines with lower normalization show 
the same congelation at z = 0.1 and z = 1 (Elbaz et al. 2007) and in the 
Millennium simulations at z = 2 (dotted line). Figure adapted from Daddi et 
al. (2007), but both SFRs and stellar masses were rescaled using a Chabrier 
IMF 



required them to have a radio continuum detection: the typi- 
cal 5ct limits of 20/iJy in our 1.4 GHz VLA maps (Morrison 
et al. 2010) correspond to SFRs of 66 Mr., yr"' at z = 1 .4 and 
SFR = 156 Mq yr-i at z = 2.0. 

We note that the six galaxies in our sample are all at 1 .4 < 
z < 1.6; i.e., in the lower redshift range probed by the BzK 
selection that extends over 1.4 < z < 2.5. This is due to the 
requirement of a radio detection coupled with the choice to 
observe the 2-1 transition of CO: such a transition has a rest- 
frame frequency of 230.538 GHz and can be observed with 
the PdBl only up to z = 1 .87. We have chosen to observe the 
2-1 transition rather than higher transitions in order to limit 
the uncertainties and corrections due to the possible presence 
of low excitation molecular gas. Based on the results by Dan- 
nerbauer et al. (2009), we would have underestimated the to- 
tal CO luminosity by a factor of ^ 2.5 if we had observed the 
CO[3-2] transition and assumed constant brightness tempera- 
ture in the CO emission. Observing CO[3-2] would have also 
resulted in substantially lower S/N ratio for the galaxies in a 
fixed integration time, a finding consistent with the non detec- 
tion of similar galaxies in CO[3-2] by Tacconi et al. (2008) 
and Hatsukade et al. (2009). 

Finally, we investigate if the requirement of the availabil- 
ity of a spectroscopic redshift prior to CO observations could 
have biased our sample. Given that most of the spectro- 
scopic redshifts in GOODS-N were derived from optical spec- 
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Fig. 9. — The reddening E{B — V) estimated from the UV slope is plot- 
ted versus the stellar mass. The large filled circles are the six CO detected 
galaxies, while smaller symbols refer to the sample of near-IR selected BzK 
galaxies (only star forming BzK galaxies, or sBzKs, are shown). Objects 
above the horizontal dashed line are not selected with the BM/BX technique 
(Adelberger et al. 2004; Steidel et al. 2004). 



troscopy, it is possible that this requirement might have biased 
our sample against optically faint galaxies or equivalently, 
given the near-lR selection of our sample, to objects with rela- 
tively low reddening by dust. Fig.|9]shows the amount of red- 
dening for the BzK galaxies, estimated from the slope of their 
UV continuum, plotted as a function of stellar mass. There is 
a clear trend between E{B — V) and stellar mass in our sam- 
ple, as reported also by Greggio et al. (2009). More massive 
galaxies are more strongly reddened, most likely a result of 
the mass-metallicity relation (Erb et al. 2006; Hayashi et al. 
2009; Onodera et al. 2010). Incidentally, we note that the vast 
majority of galaxies with M }^ 1O"M0 (including the six CO 
detected galaxies) are too red to be selected as BM/BX galax- 
ies from the UV (Fig.|9l), a result consistent with what is found 
by van Dokkum et al. (2006) at slightly higher redshifts. From 
Fig- m it is clear that the six CO-detected galaxies lie on the 
average trend of E{B — V) versus stellar mass, i.e. they do not 
show particularly bluer colors than galaxies of similar stellar 
mass. 

To summarize the characterization of our sample, we con- 
clude that we have observed a subsample of BzK galaxies 
with fairly high stellar masses and toward the lower redshift 
range of the BzK selection but with overall physical properties 
(SFRs, stellar mass and reddening) that are consistent with 
those of the general population of near-lR selected galaxies at 
z-2. 

3.2. Total IR luminosities and star formation rates 

In order to accurately derive the SFRs and IR luminosities 
of our sample galaxies we have used all the available SFR in- 
dicators. This includes the dust corrected UV luminosities, 
mid-lR continuum luminosities from 24/im Spitzer imaging 
and the VLA 1 .4 GHz radio continuum. These datasets and 
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Fig. 10. — The IR luminositie.s (Lir) for our galaxies were derived using three independent indicators: radio, UV and mid-IR. We here compare the average of 
the three indicators for each source with the individual derivations for the radio (left panel), mid-IR (center panel) and UV (right panel). In each plot we show 
a 0.17 dex en'or bar, corresponding to the scatter of all individual measurements with respect to their average. Ljr is a measure of the ongoing SFR in galaxies: 
we adopt Lir/L,., = 5f R/[Mr.-yr"'] x lO'", appropriate for a Chabrier (2003) IMF. From left to right in each panel, galaxy IDs are BzK-25536, 4171, 17999, 
16000, 21000, 12591. 



the techniques that we use to convert observed luminosities 
into SFRs (and equivalently Ljr; we adopt here Lir/Lq = 
SFR/[MQyx-^] X 10^°, appropriate for a Chabrier 2003 IMF) 
are described in detail in Daddi et al. (2007) and are largely 
based on the Chary & Elbaz template libraries (but we here 
use an improved 1.4 GHz map, Morrison et al. 2010). A com- 
parison between these (and other) SFR indicators within the 
GOODS sample of BzK selected galaxies for both individu- 
ally detected sources and through stacked sources is discussed 
in detail in Daddi et al. (2007) and shows that, on average, 
there is very good agreement between the different indicators. 

We find that the SFRs estimated using the different tech- 
niques agree generally very well also for our CO-detected 
galaxies, typically to within much better than a factor of two 
(Fig. [Tol l. In particular, no galaxy in our sample would be 
classified as a mid-IR excess galaxy according to the classi- 
fication of Daddi et al. (2007b), which could have suggested 
the presence of an obscured AGN. The largest discrepancy 
found is for BzK- 12591 where we find a radio estimated IR 
luminosity that is twice as high as both the UV and mid-IR es- 
timates. This might hint at a possible presence of an AGN, as 
also suggested by the possible [NeV] emission in the optical 
spectrum. 

We decided to adopt the average of the three estimators for 
the best measurement of Lir and, subsequently, SFR, with the 
exception of BzK- 12591 where we used the average of the 
UV and mid-IR estimates only (if we had included also the 
radio its SFR would have increased by 0.17 dex, leaving the 
remainder of our analysis basically unchanged). The uncer- 
tainty of this measure is derived using the dispersion of the 
differences between the average and each individual determi- 
nation. From this we find that our SFR should be accurate to 
~0. 17 dex. The SFRs derived from the UV (as used in Fig.|8]l 
agree with our best SFRs within 0.02 dex, with a scatter of 
0.14 dex. 

Inspection of the deep 70/im imaging from the Far Infrared 
Deep Extragalactic Survey (FIDEL, Frayer et al. 2009) shows 
that BzK-21000 is detected with a flux of 4.9 ±0.7^Jy, from 
the catalog of Magnelli et al. (2009), as noted also by Casey 
et al. (2009). The measured 70/im flux corresponds to Ljr = 
2.1 ±0.3 X 10^%) using the Chary & Elbaz (2001) library 



that we adopt, in very good agreement with our estimate in 
Tab. |4] This is an independent confirmation of our Lir and 
SFR estimates. 

Overall, the galaxies in our CO sample are actively star 
forming, with SFRs ranging from about 6QMr., yr^' for BzK- 
25536 to 4OOM0 yr^^ for BzK-12591 (Tab.l^l. 

3.3. Stellar masses from synthetic template fitting 

Multicolor photometry in the rest frame UV, optical and 
near-lR bands is available for our CO detected galaxies in 
GOODS-N, as presented by Daddi et al. (2007). The UB- 
VRlzJHK photometry, plus the photometry in the first three 
IRAC bands at 3.6/im, 4.5/im and 5.8/im is shown for all tar- 
get galaxies in Fig.fTTf". 

We used synthetic templates from the library of Maraston 
(2005) to derive stellar population properties and in particular 
stellar masses. These templates have a consistent treatment 
of the emission of AGB stars that is important for the cor- 
rect derivation of stellar masses in the distant Universe (see, 
e.g., Maraston et al. 2006). Given that all of the galaxies in 
our sample are actively star forming, we have used templates 
with a constant SFR, a large range of ages, Chabrier (2003) 
IMF and a range of metallicities (half solar, solar and twice 
solar, appropriate for the expected range in stellar masses). 
We allow for reddening using the Calzetti et al. (2000) ex- 
tinction law. We explicitly account for the mass loss due 
to stellar evolution, so that the stellar masses that we derive 
are lower by some 10-20% (depending on age; see Fig. 3 of 
Maraston et al. 2006) than the time integral of the SFR. The 
results are summarized in Tab. |4] The stellar masses range 
from3.3 x lO^Vg for BzK-25536 to 1.1 x lO'^Mg forBzK- 
25536. The fitted stellar masses are only 0.06 dex lower on 
average than those from the D04 empirical formula, with a 
scatter of 0.08 dex. This implies that Fig.[8]is relevant also for 
the stellar masses. 

We estimated the uncertainties in the derivation of the stel- 
lar masses on the basis of variations, following Avni ( 1 976) 
for the case of one interesting parameter, marginalizing thus 

Following Maraston et al. (2006) and other work we do not use the 
IRAC 8.0/im band as it appears to be often contaminated by dust emission in 
the BzK galaxies, see also Daddi et al. (2007b). 
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Fig. 1 1 . — Spectral energy distribution of our target sources showing best fitting synthetic templates from the M05 library. We consider constant star formation 
rate models with 0.5, 1 and 2 solar metallicities and allow for reddening using a Calzetti et al. (2000) law. EiTors bars in the plots have been increased to formally 

havexLw = 1- 



on the template age, reddening and metallicity. We find that, 
given our accurate multiwavelength SEDs, the formal (statis- 
tics) la error on the stellar masses are of order of 0.10- 
0.15 dex. 



4. STAR FORMATION EFFICIENCIES 



From the observed CO[2-l] fluxes we have de- 
rived CO luminosities L^q spanning a range of 1.2- 
3.2x10'^ K km s"' pc^ (Tab. H see Solomon & van 
den Bout 2005 for the conversion formula in the adopted 
units). We here have applied a 16% correction to convert 
luminosities derived from CO[2-l] to Lqq (defined relative 
to the fundamental CO[1-0] transition) based on the BzK 
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Fig. 12. — CO luminosities (L'qq) versus the bolometric IR luminosities (Lir, left panel) and versus the stellar masses (right panel). The dotted line in each 
panel is the best linear fit. The average ratio and dispersion is reported in the bottom of each panel (units are K km s^' pc^ Lp' for the ratio in the left panel and 
K km s"' pc - M Tt' in the right panel). The red solid line in the left panel correspond to the fit to BzK galaxies at z = 1.5 and spiral galaxies, as shown in the left 
panel of Fig.[T3] In both panels galaxy IDs are, from top to bottom: BzK-12591, 21000, 4171, 17999, 16000, 25536. 



Table 4 

Physical properties of the CO detected galaxies 



Source 






''CO 


SFE 








Mdy„« 


aco 


gas 


t 8 

'rot 


fSFR 




Mr; yr-' 


10'2l.: 


lO'^Kkms- 


' pc^ L,;;K"'km-' spc-- 


1O"'M0 


lO'^M,; 




io'"m,3 


MoK^'km"' s pc" 


- Gyr 


Gyr 


Gyr 


BzK-4171 


103 


1.0 


2.1 


49 


4.0 


7.7 


0.66 


8.2 


3.9±1.2 


0.74 


0.09 


0.39 


BzK-21000 


220 


2.2 


2.2 


98 


7.8 


8.1 


0.51 


9.3 


2.7±1.4 


0.37 


0.17 


0.36 


BzK- 16000 


152 


1.5 


1.6 


93 


4.3 


5.9 


0.58 






0.35 


0.34 


0.28 


BzK-17999 


148 


1.5 


1.7 


85 


3.9 


6.3 


0.62 


7.3 


4.0±1.3 


0.43 


0.13 


0.27 


BzK- 12591 


400 


4.0 


3.2 


124 


11 


12 


0.52 






0.29 


0.09 


0.27 


BzK-25536 


62 


0.6 


1.2 


53 


3.3 


4.2 


0.56 






0.68 


0.21 


0.53 



Notes: 

1 : SFRbest is the average of dust corrected UV, mid-lR and radio deiived SFRs. Typical en'ors are at the level of 0. 17 dex (see text). 
2: we use Ljr/L,, = 5f R/[M.,yr"'] x 10'°, adapted from Kennicutt et al. (1998) using a Chabrier IMF. 
3: we apply a coirection of 16% to account for the CO[2-l] transition being slightly sub-thermally excited; Dannerbauer et al. (2009). 
4: computed from L'qq, assuming oco = 3.6 as deduced from our observations 
5 : /ga.s = ^gas / (Wgas + Asians) ; the average error is 9% (including also the uncertainty on aco 
6: estimated for galaxies resolved in CO emission, except BzK-16000 that is seen nearly face-on; errors on Mdy„ range from 30 to 50%. 

7:%as=Mg,,/SFR 
8: /rot = 27rrj/ (i'fwhm/2) is the rotation time at the half light radius 
9: rsFR=M„a,/SFR 



excitation results in Dannerbauer et al. (2009). Given the flux 
measurement errors and the uncertainty in this conversion, 
we estimate typical errors in Lqq to be at the level of 0. 1 dex. 

We compare the CO luminosities to the IR luminosities and 
stellar masses of our sample galaxies in Fig. [12] We find that 
the CO emission correlates well with both Lir (hence SFR) 
and stellar mass (we exclude at the 2.5a [3a] confidence level 
that the CO emission does not correlate with SFR [mass]). 
The ratio of L[^q/Lir (i.e., SFE) has a scatter of only 0. 16 dex, 
while the individual values are all within a factor of 2.5. The 
scatter is in fact similar to the error of Ljr only, which possi- 
bly dominates the total error (this would suggest even lower 
intrinsic scatter in Ljr/ Lqq). The ratio of L^o/^stars has an 
even lower scatter of 0.10 dex. This is driven by the fact that 
the formal errors on stellar masses are lower than those on 
IR luminosities. This correlation of the CO luminosity with 
both the SFR and stellar mass is not surprising, given that 
the six galaxies lie on top of the stellar mass-SFR correlation 



and their SSFR has a scatter of only 0.12 dex. However, it 
is remarkable that we can now probe for the first time that 
the narrow range of SSFR corresponds to a similarly narrow 
range of SFEs, i.e. is directly related to the gas properties. 

In Fig. [13] we compare the CO and IR luminosities of the 
six CO-detected BzK galaxies with other cosmologically rel- 
evant galaxy populations: SMGs (Greve et al. 2005; to which 
we add more recent results of Daddi et al. 2009ab and Frayer 
et al. 2008), quasars (QSOs; Riechers et al. 2006 Solomon 
& van den Bout 2005), local ULlRGs (Solomon et al. 1997) 
and local spirals taken from the HERACLES survey (Leroy et 
al. 2008; 2009) and from the Virgo Cluster survey of Wilson 
et al. (2009). We notice that in Daddi et al. (2008) we used 
the results of Yao et al. (2003) to define a bona fide compar- 
ison sample of local spirals, following Solomon & van den 
Bout (2005) and Greve et al. (2005). While the results and 
conclusions of this paper would remain the same if we were 
still using the Yao et al. (2003) sample, we here prefer the 
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Fig. 13. — Comparison of the CO and IR luminosities of tlie CO-detected BzK galaxies (red filled circles) to various local and distant samples: SMGs (Greve 
et al. 2005 and Daddi et al. 2009: blue empty circles; Frayer et al. 2008: blue empty square), QSOs (Riechers et al. 2006; Solomon & van den Bout 2005; green 
empty triangles), local ULIRGs (Solomon et al. 1997: crosses), local spirals (Leroy et al. 2008; 2009: filled squares; Wilson et al. 2009: filled triangles). For 
local spirals we associate a Hubble-flow redshift (right panel) based on the actual distance. The solid line in the left panel shows the best fitting L'qq versus Lir 
relation to the combined sample of local spirals and distant BzK galaxies. The dotted line shows the same relation with normalization higher by 0.5 dex. 



more recent samples as the Yao et al. sample was selected in 
the submm. Also, we here chose not to adopt the IR lumi- 
nosities for SMGs estimated by Greve et al. (2005) but use 
a lower conversion factor between 850/xm flux and IR lumi- 
nosity to account for the evidence that their SEDs are colder 
than previously thought (e.g.. Pope et al. 2006; 2008; Kovacs 
et al. 2006; Daddi et al. 2009). In addition, we exclude lower 
luminosity objects from the Solomon et al. (1997) sample re- 
quiring LjR > 0.7 X 10^^ to obtain a sample of ultra-luminous 
IR galaxies with IR luminosities similar to the BzK galaxies. 

Based on our sample of six galaxies, we thus confirm and 
extend the result of Daddi et al. (2008) that the typical SEE 
of massive BzK galaxies is similar to that of spirals and sig- 
nificantly lower than that of SMGs. We find that the aver- 
age SEE (in units of LqK^^ km^^ s pc^^; the errors given 
in the following are the uncertainties of the mean) of BzK 
galaxies is 84 ± 12, that of SMGs is 390 ± 60, for QSOs we 
find 590 ± 100, for local ULIRGs we have 250 ± 30 and local 
spirals have 48 ± 7. Therefore, the typical BzK galaxy has 
about 80% higher SEE of that of a local spiral, but individ- 
ual values are well within the range measured for local spirals 
even though BzKs are a hundred times more luminous. Local 
ULIRGs have SEEs that are three times larger on average than 
BzK galaxies. Although there is some overlap in the spread 
of individual targets, one can see that the two distributions do 
not really overlap in the Lir versus L^q plane (Eig.[T3]l. 

Eor SMGs/QSOs the average SEE is 4.5/7 times larger than 
that of an average BzK galaxy. Although at least one case 
is known of an SMG with a SEE at the level of that of BzK 
galaxies (Erayer et al. 2008), for the vast majority of the sam- 
ple the distribution of values are fairly disjoint. We note that, 
in contrary to our BzK galaxies, the CO luminosities of SMGs 
are derived from high-7 transitions, and could be underesti- 
mated if a substantial amount of cold gas is present. Such a 
correction could indeed be important in some cases (Greve et 
al. 2003; Carilli et al. 2010) but the difference still appears to 
be quite significant if a correction was applied using the av- 



erage CO SEDs of SMGs (Weiss et al. 2005; 2007) or local 
ULIRGs (Papadopoulos et al. 2007). 

Given the comparable SEE of local spirals and BzK galax- 
ies, and the evidence that indeed the BzK galaxies are gas rich 
disks in the distant Universe, we can formally consider them 
as a homogeneous population and derive the best fitting Lqq 
versus Ljr relation for the combined population. We find the 
relation: 

log(LiR) = 1.13 x(logL^o) +0.53 (1) 

with Lqq in units of K km s^ ' pc^ and Lir in units of Lq. Both 
parameters in the fit are determined with an accuracy of 0.07. 
The relation appears to fit both the sample of BzK galaxies 
(see also Eig. fT2l i and local spirals well, with a root mean 
square scatter of 0.24 dex. The local ULIRGs and SMGs are 
systematically offset from the fit, but are consistent with a 
relation having a similar slope but with a factor of 3-4 higher 
normalization. 

Similarly, Eig.[T3]shows that our BzK galaxies should have 
been factors of 3-4- fainter in their CO fluxes had they fol- 
lowed the trends previously established for equivalently lumi- 
nous galaxies such as local ULIRGs and SMGs, as discussed 
in the introduction. Had that been the case, it would have 
taken ^10 times longer integrations to detect these sources, 
which would have made detections impossible until ALMA is 
in full operation. 

5. GAS KINEMATICS IN CLUMPY DISKS: NUMERICAL 
SIMULATIONS 

A critical piece of information that our observations can 
provide is the derivation of the dynamical mass of the sys- 
tems through the observed gas kinematics and spatial extent. 
In today's disk galaxies, the cold gas component is slowly ro- 
tating and directly traces the dynamical mass, which relates 
to the circular velocity squared and the radius of the system 
(e.g., de Blok et al. 2008). However, gas motions in high- 
redshift disk galaxies are more complex. In particular, larger 
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Fig. 14. — A few examples of our simulations of clumpy disk galaxies, 
based on the results of Boumaud, Elmegreen & Elmegreen (2007). Top pan- 
els show the distribution of stellar light as observed in the UV (tracing mainly 
star formation), bottom panels that of the molecular gas. The distribution of 
star formation and optical morphologies are dominated by giant clumps, and 
our models include a variety of dumpiness and bulge fractions consistent 
with observed star-forming disks at high redshift. These snapshot show a 
40 X 40 kpc area (5" x 5" at z = 1 .5). 
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Fig. 15. — Calibration of the relation between the RxV^ estimator cor- 
rected for inclination effects (y-axis) and the effective amount of mass con- 
tained within the half light radius (;-<,, x-axis), based on clumpy disk models. 
Sequences of points at fixed dynamical mass correspond to the same model 
observed with different inclinations and azimuths. The solid and dashed fines 
sfiow tfie best fitting relation (Eq. |2j and the Ict scatter on the relation, re- 
spectively. 



molecular gas reservoirs are present and optical morphologies 
are dominated by giant star-forming clumps. This is seen in 
our sample (Fig.|5]) as well as in the majority of star-forming 
disks at high redshift (Elmegreen et al. 2007). These gi- 
ant clumps are up to ^ lOOOx more massive than the star- 
forming clouds of today's spirals (Elmegreen et al. 2009b). 
Such gas-rich clumpy disks can have significant non-circular 
motions (Immeli et al. 2004; Bournaud et al. 2008) and show 
strong turbulence in their ionized gas (e.g., Genzel et al. 2008; 
Forster-Schreiber et al. 2009). It is likely that such motions 
are also present in their cold gas (CO) component (Bournaud 
& Elmegreen 2009; Dekel et al. 2009b). dumpiness and tur- 
bulence influence the kinematics, sizes, velocity fields, and 
dynamical mass estimates of distant galaxies (Bournaud et al. 
2008;Burkertet al. 2009). 

To recover dynamical masses accurately from CO velocities 
and to estimate associated uncertainties, we used numerical 
simulations of gas-rich clumpy galaxies based on the work of 
(Bournaud, Elmegreen & Elmegreen 2007; hereafter BEE07). 
These simulations start as gas-rich, gravitationally unstable 
disks that rapidly evolve into clumpy disks with a phase space 
distribution similar to cosmological simulations of clumpy 
disk formation (Agertz et al. 2009; Ceverino et al. 2009) 
and consistent with observed clumpy galaxies (see BEE07 
and Bournaud et al. 2008). The initial conditions (f — 0) of 
the six new models used here were scaled to match the typical 
masses and sizes of the observed galaxy sample, and are as 
follows: stellar masses ofexlO'^Mg, initial half mass radii 
of 5 kpc, gas fractions /g = (Mgas/Mgas+Mstars) of 30, 50, and 
70%, and dark matter fractions /dm = (A^daik / (A/dark +A^gas + 
A^stars) of 15 and 35% within the effective radius. The out- 
put of the simulations were retained after 250, 450, 650Myr, 
which, in the evolutionary sequence studied in BEE07, cov- 
ers all stages from very clumpy bulgeless morphologies to 
smoother bulgy systems, consistent with our observed sam- 
ple (Fig. [141 compare to Fig.|5]). In each case the system was 
observed with 30" and 60" inclinations and at three different 
azimuths, giving a total of 108 "independent" projections. 

Taking advantage of the fact that the actual mass encom- 



passed at the effective radius can be directly measured in 
the simulations, we compare its value to a crude estimate of 
the dynamical mass, (vfwhm/2)^ x r^/ (G sin^ /), in order to 
derive a robust method to estimate dynamical masses from 
high-redshift CO spectra. For each projection of the simu- 
lations, we derive the spatial distribution of "molecular" gas 
component. To do so, we only consider gas in regions where 
the density is higher than 100 atoms cm^^ We measure 
the effective radius not as the true half-light radius but as 
the effective radius of a Sersic fit to the radial profile distri- 
bution, just as in our observations. We derive synthetic CO 
spectra from which we estimate the FWHM velocity vfwhm- 
We also produce UV stellar images from the simulations to re- 
produce deprojection uncertainties and compare UV and CO 
sizes (see below). In this way, the (vfwhm/2)^ x r^./ (G sin^ /) 
ratio is computed with the same definitions, assumptions and 
possibly biases as in the observations. The main result of this 
exercise is shown in Fig. [15] We find the following best fitting 
relation to the actual mass inside re'. 

re X fvFWHM/2)^ 

M{r <re)^l3x- — ^ ^^^f-' ' ±12.5% (2) 
G sin / 

The correction factor of 1.3 includes various effects: 

• the flattening of the baryonic mass distribution in a disk 
(Binney & Tremaine 2008). This has a small effect 
(around —5%) on the total (+30%) correction factor. 

• a possible bias on estimating the projection angle / 
because of the dumpiness of these galaxies. We de- 
projected simulated UV images just like observations, 
so that any uncertainty and bias on / are included in 
our results (Eq. |2]and the scatter of Fig. [TSl l. We actu- 
ally found that the inclination / retrieved from the UV 

' ' This selection actually includes almost all the gas mass in the models, 
since the gas fraction is high and the density is high almost everywhere. 
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Fig. 16. — Four CO spectra extracted from our simulations. The grey 
histogram shows the initial (f = 0) uniform disk, the black lines and yel- 
low histograms correspond to evolved stages in the simulations (i.e., realistic 
clumpy disks). Because of the discrete sampling of the velocity field related 
to the dumpiness of the gas, and the individual motions of each clump, the 
predicted spectra are more irregular than for z = spiral galaxies. Overall, 
the typical double-horn profile is preserved, but there can be a strong central 
dip at zero velocity (a,c, resembling the observed spectrum of BzK-4171 in 
Fig. |3) or, on the other hand, relatively flat spectra with little zero velocity 
cenfi'al depression (b, as in BzK-17999). Asymmetrical spectra (a,d) are also 
found, as observed in BzK-12591 and BzK-25536. The velocity width of 
these spectra is smaller than for the initial smooth rotating disk. This corre- 
sponds to clump evolution and migration with rotation velocities below the 
circular velocity, and the associated turbulent pressure support that compen- 
sates for the lack of rotation. This lies at the origin of the calibration factor in 
Eq.|2] 



images tends to slightly overestimate'^ the actual pro- 
jection angle (which is known in simulations). This is 
also a small effect, around +10% on M{r < re). 

• The fact that gas motions largely differ from circular 
orbits is the main effect, around +25% on M{r < re). 
The giant and massive gas clumps interact with each 
other, undergo dynamical friction and migrate radially 
(BEE07, Dekel et al. 2009b). The presence of a high 
velocity dispersion a decreases the observed rotation 
velocity Vrot compared to the actual circular velocity 
Vcirc tracing the potential, with v^jj.^, ~ v^^j + 3ct^, where 
(7 is the one-dimensional gas velocity dispersion (ob- 
served line-of-sight dispersion) and the factor 3 is for 
isotropic gas turbulence in three dimensions. The effect 
is also known as 'asymmetric drift' in stellar kinemat- 
ics, and the equivalent of a turbulent pressure support in 
addition to the rotational support for a gas disk (Burkert 
et al. 2009; Elmegreen & Burkert 2009). The velocity 
dispersion a of the CO component cannot be directly 
measured given the current observations, but should be 
several tens of km s^' based on Ha data (e.g., Forster- 
Schreiber et al. 2009) and observed clump sizes (which 
is set by the Jeans mass, Bournaud & Elmegreen 2009). 
The theoretical correction factor for the mass expected 
from a disk with a circular velocity of about 200 km s^' 
is then consistent with the +25% correction from our 
simulations. 

Our simulations also have the advantage of including the 

For instance, a face-on disk (;' = 0) could have a large fraction of its 
light distributed in a few main clumps, leading to an axis ratio b/a < I and 
hence a derived inclination > 0. 



dumpiness of the gas resulting from an irregular sampling of 
the underlying velocity field and a dispersion in the dynami- 
cal mass measurements - a given model has different velocity 
widths vpwHM when we observe it from different azimuths. 
The resulting uncertainty on the estimated mass is neverthe- 
less relatively small, about 12.5% (this uncertainty accounts 
for the method used and dumpiness of disks). In summary, 
the presence of giant clumps does not significantly hamper 
the possibility to estimate accurate dynamical masses. 

The models were started with similar gas and stellar sizes, 
and we checked that the half light radii derived from the UV 
stellar emission and from the CO emission remain similar dur- 
ing the evolution (Fig.O. This is in agreement with what is 
found for our sources (see Section |2]l as well as local disks 
(Leroy et al. 2008). This implies that, referring to the dynam- 
ical mass within the half light radius, the quantity computed 
according to Eq. |2] should be compared to half of the total 
amount of stellar mass and gas mass in the system, plus the 
amount of dark matter eventually present within r^, i.e.: 

M{r < re) = 0.5 X (M^tar +Mgas) +Mdark(r < r,) (3) 

Applying this result to observations, we derive M{r < re) 
from Eq. |2] subtract the known stellar mass Mjtar.tot and an 
estimate of M(iaiii('" < r,.) to finally derive the gas mass Mgas.tot 
and compare to the observed CO fluxes. 

Note that we calibrated the measurements (Eq. |2] and 
Fig.fTsll at a radius r = re and using the vfwhm velocity, in or- 
der to match observational measurements. Dynamical masses 
estimates at different radii should be similar, as the main con- 
tribution to these factors in Eq.|2] namely dumpiness and high 
turbulence, are present all over high-redshift disks. 

The numerical simulations spanned gas fractions from 30 
to 70%. Reproducing clumpy structures similar to those ob- 
served (Fig. |5] and |2]i requires high gas fractions and disk 
masses (Bournaud & Elmegreen 2009), and we will show 
that our observations support a high molecular gas fraction 
in these galaxies. In any case, kpc-sized clumps are observed, 
which imply a large Jeans length and a high turbulent speed 
(velocity dispersion). As the underlying physics in Equation|2] 
is based on the gas velocity dispersion, it should remain valid 
independent of the actual gas fraction. 

5.1. Explaining the shapes of the observed CO spectra 

An application of the numerical models is to produce syn- 
thetic CO spectra and compare their shapes to the observed 
ones. Indeed, the observed spectra have some particularities 
that would not be expected for regularly rotating disks such as 
seen in z = spirals. 

For example, the CO spectrum of BzK-4171 (Fig. O has 
a very marked and quite symmetric double horn profile with 
a pronounced central depression in the center, falling to only 
15% of the peak flux. Gas spectra in local galaxies show dou- 
ble horn profiles but the emission at zero velocity never falls 
below 50% of the peaks (e.g., Walter et al. 2008). The ob- 
served feature in BzK-4171 can not be easily explained by 
assuming that the molecular gas is in a ring: removing molec- 
ular gas from the system barycenter would still leave substan- 
tial amount of emission at zero velocity, coming from regions 
where the gas rotation velocity is perpendicular to the line 
of sight. For instance, even a pure thin ring rotating with 
V/ct = 10 would have a zero velocity emission at 42% of the 
peaks. Fig. |3] shows that also BzK-12591 and BzK-25536 
have significantly asymmetric spectra, with a flux in one horn 
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Fig. 17. — Constraints on the CO luminosity to gas mass conversion factor (aco). The analysis was possible for three galaxies for which we plot gas masses 
versus CO luminosities (left panel) and their ratio (ctco) versus the dynamical masses (right panel). The dashed lines show the la range on aco derived by 
combining the measurements for the three galaxies. The dotted lines show the conversion factor for the Milky Way galaxy (top) and for local ULIRGs (bottom). 
From left to right galaxies are: BzK-17999, 4171, 21000. 



much higher than in the other horn. The spectrum of BzK- 
17999 is almost flat around the central velocity. 

The synthetic CO spectra from our numerical simulations 
shed light on this situation. Spectra qualitatively similar to 
BzK-4171 are found, as well as spectra where the zero ve- 
locity level is quite high, and double horn profiles that are 
strongly asymmetric (Fig. \T6[ . The explanation is that much 
of the CO emission comes from a few big clumps, and there- 
fore is discretized. At particular projections, there can be 
a lack (or an excess) of zero-velocity emission, when few 
(many) clumps have velocities perpendicular to the line of 
sight. Significant asymmetries arise when the clumps are not 
uniformly distributed between the blueshifted and redshifted 
parts of the disk. 

The successful reproduction of the observed spectra give 
confidence that our models are correct, at least to first or- 
der. More generally, both observations and modeling show 
that turbulent and clumpy disks at high redshift can have inte- 
grated spectra that are distinct from what is seen in the local 
Universe. 

An interesting characteristic in the simulated spectra is that, 
in spite of their disturbed shapes, the velocity spread remains 
about the same in the two horns (red- and blue-shifted sides). 
One of the two sides can have a much higher flux, but not a 
much larger velocity spread. This is seen also in the observed 
spectra. 

6. ESTIMATING GAS MASSES IN NORMAL HIGH REDSHIFT 
MASSIVE GALAXIES 

Based on the results shown in Fig.[T5]and using Eqs.|2]and 
[3] we proceed to estimate the dynamical mass within of 
the CO-detected BzK galaxies, based on the vfwhm velocity 
from CO[2-l] spectra, the stellar half light radius (that corre- 
sponds to the CO half light radius), and the inclination / de- 
rived from the b/a axis ratio seen in the ACS imaging. Values 
and uncertainties are given in Tab.|2]and[3] We limit our anal- 
ysis to the CO-detected galaxies that were resolved in CO. We 
did not estimate a dynamical mass for BzK- 1 6000 because the 



correction due to its inclination is quite uncertain {bja ~ 1). 

In order to estimate the gas mass following Eq. |3] we ac- 
count for a dark matter fraction within of 25% (see also 
next section). Disk galaxies are known to be dominated, at 
their effective radius, by baryons more than by dark matter 
(Bosma 1981; de Blok et al. 2008; Trachternach et al. 2008 
- dark matter dominates rotation curves only at much larger 
radii). This implies that the dark matter fraction is well below 
50%, and there is no expectation for a major change with red- 
shift (see also Forster-Schreiber et al. 2009). The structural 
properties of z 2 galaxies also indicate that, within their vis- 
ible, star-forming extent, they have baryonic masses and disk 
masses that should account for at least two thirds of the total 
masses (e.g., Bournaud & Elmegreen 2009). 

We find dynamical masses within for the three galaxies 
in the range of 7-9x IO'^Mq with typical uncertainties at the 
level of 30%. Subtracting the stellar and dark mass from the 
dynamical mass (Eq. [3]l, we derive a total gas mass of order 
6. 1-8.3 X lO^^M© for the three galaxies. After propagating all 
the estimated uncertainties in our measurements and adding 
in quadrature the 12.5% uncertainty associated to Eq. |2l the 
uncertainty in the total gas mass estimate is about 30% for 
BzK-4171 and 17999 and 50% for BzK-21000. 

These gas mass estimates are derived independently from 
the observed CO luminosities. Thus taking the ratio of the gas 
mass estimates and the CO luminosities, we can derive a con- 
version factor aco = Afgas/i'co f^'" ^^^^ system (see Tab. |4] 
and Fig. 117b. and find a combined result of 

aco = 3.6 ± 0.8 M0(K km s^'pc^)-' (4) 
Although the uncertainty in the individual aco esti- 
mates is significant, their aco value are well above the 
0.8 M,;j (K km s^' pc^)^' factor that is generally applied to 
ULIRGs and SMGs (Downes & Solomon 1998; Solomon & 
van den Bout 2005; Tacconi et al. 2006; 2008) with individual 
values for ULIRGs typically in the range of 0.3-1.3. The in- 
dividual values are also somewhat lower than the MW factor 
of 4.6 Mq (K km s^^ pc^)^', although still consistent within 
^1(7 for each object. 
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Fig. 18. — Gas fractions for the six CO detected BzK galaxies, derived using a MW-Iike aco = 3.6 ai'e compared to stellar masses (left panel) and SFRs 
(right panel). The solid and dotted horizontal lines show the average value and the la scatter. The red curves show the expected trend assuming at face value 
the correlation between L^q and Ljr (hence SFR) and the stellar mass-SFR correlation at z 1.4 (see discussion in the text). The black filled squares in the left 
panel show the molecular gas fraction for local spirals with M > 2x lO'^A/g, from Leroy et al. (2008), including H2 and helium (for our BzK galaxies we expect 
that most of the hydrogen is molecular). From bottom to top the galaxy IDs are: BzK-21000, 12591, 25536, 16000, 17999, 4171. 



Combining the three galaxies allows us to exclude the 
ULlRG-like conversion factor at the 3.7(7 level, while the 
MW-like value is within 1. 3a of our estimate. 

This conversion factor in turn allows us to estimate the gas 
content in the BzK galaxies for which a direct estimate was 
unfeasible due to larger observational uncertainties. The gas 
masses and fractions in these three galaxies turn out to be sim- 
ilarly high as in the three galaxies with dynamical estimates 
(Tab. Si. 

We notice that our aco traces the total amount of gas, in- 
cluding molecular and atomic hydrogen as well as helium 
in their respective proportions within the galaxies half light 
radii. However, we expect the vast majority of the hydrogen 
to be molecular given that the high observed densities and 
expected pressures in the interstellar medium of these galax- 
ies (see, e.g.. Blitz & Rosolowsky 2006; Leroy et al. 2008; 
Obreschkow & Rawlings 2009). This is also consistent with 
hydrodynamic simulations where more than 90% of the gas 
mass is molecular, at densities above 100 cm^^ and tempera- 
tures below ^ 100 K (Delaye, Bournaud et al., in preparation). 

In summary, our analysis shows the existence of large 
molecular gas masses in our targets, and a high CO-to- 
gas mass conversion factor, similar to the one found in the 
Milky Way. We can estimate the gas mass fractions /gas ~ 
Mgas/{Mgas +M stars)- For consistency, we apply the average 
conversion factor of aco = 3.6 to all observed CO luminosi- 
ties. The gas fractions range from 51 to 66%, with the aver- 
age being 57% with a dispersion within the sample of of 6%. 
These galaxies clearly have very high gas fractions, with gas 
accounting for masses that are similar or greater than those of 
the stars. 

The gas fraction appears quite constant within the sample 
(Fig. [TST l. Given the correlation between L'qq (and thus gas 
mass) and Lir (or SFR) with a slope above unity, and the cor- 
relation between stellar mass and SFR with a slope slightly 
below unity, one could expect a gas mass fraction slowly de- 
creasing with increasing stellar mass and SFR. The data are 



consistent with such a weak trend. 

For comparison, the total gas fractions (including HI, H2 
and helium) in local spiral galaxies with the same mass is 
about 20% (Leroy et al. 2008), and about 7% if including 
only H2 and helium (i.e., restricting to the stellar disk). We 
thus observe a very strong evolution in the gas fractions of 
disks at fixed stellar mass. Assuming that most of the gas 
within re in BzK galaxies is molecular, this corresponds to an 
average increase of molecular gas content of a factor of 18, at 
fixed stellar mass, from z = to 1 .5. This is very similar to the 
increase of the SSFR (or the normalization of the stellar mass- 
SFR relation) between z = and 1.5. Hence the increase in 
SSFR can be consistently explained by the increase of molec- 
ular gas content inside star forming galaxies. 

Dividing the gas masses by the SFRs we estimate gas con- 
sumption timescales in the range of 0.3-0.8 Gyr. The SFRs 
can in principle continue for even longer timescales, once ac- 
counting for gas accretion. This confirms that star forma- 
tion can go on for much longer timescales than expected in 
rapid, generally merger driven bursts. The (rotation) dynami- 
cal timescale at the half light radius is frot — 27rre/ (vfwhm/2), 
and ranges over 0.1-0.3 Gyr for our galaxies. On average 
the gas consumption time thus corresponds to 3.5 rotation 
timescales for the disks at r^,. 

6.1. Discussion of systematic uncertainties 

It is important to investigate possible systematic effects that 
could be biasing our dynamical mass analysis and the result- 
ing derivation of the aco conversion factor Generally speak- 
ing, these CO detected BzK galaxies have bolometric IR lu- 
minosities similar to those of local ULlRGs, but we advocate 
a much larger conversion factor. Several reasons make this 
quite plausible as these galaxies behave like disk galaxies in 
many regards - see discussions in Daddi et al. (2008), Dan- 
nerbauer et al. (2009) and in the next section. However, we 
should investigate if there are issues that could have led to an 
overestimation of the conversion factor. 

A first possible issue deals with systematics in the deriva- 
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tion of the dynamical masses. Our analysis is calibrated on 
clumpy disk numerical modeling, which seems well justified 
given that our sources are observed to be highly clumpy. If 
we were to neglect the factor of 1.3 in Eq. |2] the combined 
constraint on aco would get smaller, but still 3 times higher 
than the ULlRG-like value. Nevertheless, this correction fac- 
tor is unlikely to be a numerical artifact, as it is theoretically 
justified from the high velocity dispersions observed in high- 
redshift gas disks (Burkert et al. 2009). 

If we interpreted our sources as merging galaxies rather 
than clumpy disks, the implied dynamical masses would have 
grown even higher (see, e.g., Tacconi et al. 2008), and so 
would the conversion factor. 

We might also consider systematic effects in the determina- 
tion of the stellar mass. We have used a bottom-light IMF 
that is the best choice according to observations in the lo- 
cal Universe (Chabrier 2003; see also Kroupa et al. 2002). 
Several authors have suggested that more top-heavy IMFs 
might be appropriate in the distant Universe (Dave 2008; van 
Dokkum 2008), which would further decrease the estimated 
stellar masses and hence increase the conversion factor. 

There remains the possibility that these galaxies could con- 
tain large reservoirs of old stars, formed at very high redshifts, 
that could be missed because the light is dominated by young 
stars. In practice, we would need to postulate hidden stellar 
mass reservoirs that should increase the total stellar masses 
by factors of 2-3 in order to bring our estimate of aco to the 
ULlRG-like level. However, here we use observations up to 
the near-IR rest-frame (from IRAC) that should correctly re- 
trieve the masses of old stars without such a strong bias. Also, 
these putative old stellar reservoirs would most likely be in 
bulges and halos, which actually can not contain most of the 
stellar mass around z ^ 2 disks according to the modeling of 
Bournaud & Elmegreen (2009). In addition, the existence of 
the stellar mass-SFR correlations suggests that the appropriate 
star formation histories for high redshift galaxies have SFRs 
that are rising rapidly with time (e.g., Renzini 2009) so that 
the galaxies that we observe formed most of their stars quite 
recently. 

Another assumption that we have made is a dark matter 
fraction of 25% inside the effective radius. There is of course 
a plausible range around this value. If we change the estimate 
of the dark matter fraction to 20% (30%) the average aco 
would increase (decrease) by about 10%. A survey of local 
spirals by Pizagno et al. (2005) shows a total-to-stellar mass 
ratio around two within 2.2 times the exponential disk scale- 
length for galaxies with baryonic masses ~ 10" Mf), which 
for an exponential disk corresponds to a total-to-stellar ratio 
of 1.4 within the half-light radius, i.e. a dark matter fraction 
of 20-25%. Padmanabhan et al. (2004) find similar dark mat- 
ter fractions in massive galaxies and emphasize the fact that 
higher fraction are found only in low-mass and low-density 
galaxies. The recent mass models of the Milky Way by Xue et 
al. (2008) correspond to a dark matter fraction of 23% within 
the half-light radius, using NEW profiles for the halo. There 
is no known reason for a major change with redshift'^. The 
observed objects could turn into elliptical galaxies by redshift 
zero, but the dark matter fraction still remains limited to less 
than 25%, and likely 15%, within the half-mass radius (e.g., 
Gnedin & Ostriker 2000). Thus, our choice of 25% lies on 

Adiabatic contraction of halos may actually increase the dark matter 
mass enclosed by a given radius with time, but the effect should be weak 
between z = 1 .5 and z = 0. 
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the conservative side of the plausible range. A dark matter 
fraction of 60%, which would lower our aco to the ULIRG 
value, appears to be very unlikely. 

7. IMPLICATIONS ON STAR FORMATION IN DISK GALAXIES AT 
HIGH REDSHIFT 

7.1. Massive, extended reservoirs of low-excitation 
molecular gas 

CO[2-l] detections of two near-lR selected (BzK) galaxies 
at z = 1 .5 were presented in Daddi et al. (2008). We have now 
observed a total of six galaxies in CO[2-l] and detected all of 
them. This clearly supports the ubiquity of massive molecular 
gas reservoirs in these galaxies, a result confirmed by Tacconi 
et al. (2010). 

In our observations we were able to spatially resolve ro- 
bustly the CO[2-l] emission in this kind of galaxies and derive 
typical CO emission EWHM in the range of 6-11 kpc. This 
is a factor of 2-3 larger than what is found in SMGs (Tacconi 
et al. 2006; 2008; Bouche et al. 2008). The presence of ex- 
tended gas disks is quite consistent with our previous finding 
based on the low CO excitation in these objects (Dannerbauer 
et al. 2009). These findings are indicative of colder and/or 
lower density gas than what is found in SMGs, a situation 
more similar to what is found in local spiral galaxies. 

7.2. Normal star-forming galaxies at high redshift as 
extended, clumpy rotating disks 

Our sample is based on the BzK-selection technique and 
we have shown it to be representative of "normal", typical 
star-forming galaxies at high redshift. We now summarize the 
many pieces of evidence that these CO-detected BzK galaxies 
are large disk galaxies rather than on-going violent mergers. 

Eirst, we clearly detect double peaked profiles and the ve- 
locity width in the blue- and redshifted components is very 
similar. This would be a highly contrived coincidence if this 
was due to a merger: the components would need to have 
equal mass and be observed at similar inclinations. 

Second, the UV rest-frame morphology disfavors on-going 
major mergers: in no case do we see evidence for a double- 
galaxy encounter. The irregular and clumpy morphology 
could lead one to speculate that this is due to the presence 
of multiple minor merging events, but in all likelihood we are 
looking at clumps in single massive disks (Elmegreen et al. 
2007, 2009a, Bournaud et al. 2008, Genzel et al. 2008). 

Other evidence for rotating disks include the position of 
these sources in the velocity-size plane. Bouche et al. (2007) 
used Ha resolved spectroscopy of optical and near-lR se- 
lected galaxies and identified a dichotomy between SMGs and 
extended disks in such velocity-size diagrams: SMGs have 
much higher velocities and smaller sizes. Using CO velocity 
widths and CO sizes for both BzK galaxies and SMGs, we 
confirm a dichotomy between the two classes: in the notation 
of Bouche et al., our CO-detected galaxies have 'maximum 
rotational velocities' — vfwhm /^.-liS ^ 200 km s^' and ex- 
ponential disk scale lengths r^j > 3 kpc, while the SMGs in 
Bouche et al. have v^ > 200 km s^' and < 3 kpc. The lat- 
ter are better understood in terms of merging systems, while 
the former appear consistent with being non-interacting disk 
galaxies (or interacting only with low-mass satellites). 

Further evidence that we are looking at disk galaxies comes 
from the CO properties. We have already discussed the large 
spatial sizes of the CO reservoirs, and the low gas excitation. 
Furthermore, the SEE in these galaxies is considerably lower 
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than that of both SMGs and local ULIRGs, suggesting a dif- 
ferent star formation mode. The gas depletion timescales are 
of order of 0.5 Gyr and imply a high duty-cycle (see also 
Daddi et al. 2005; 2007; 2008) - this is quite different from 
the timescales expected for a merger-induced event that would 
result in a ^ lOOMyr starburst (Mihos & Hernquist 1996; 
Greve et al. 2005; Di Matteo et al. 2008). As can be seen 
from Fig. [T3] SMGs typically also have depletion timescales 
of <100Myr. 

7.3. Origin of massive disks at high redshift 

The fact that we are presumably looking at non-interacting 
rotating disks instead of on-going mergers does not in itself 
imply that these massive disks were not assembled by past 
mergers. Remnants of major mergers can have massive rotat- 
ing gas disks, especially when the gas fraction is high (e.g., 
Robertson & Bullock 2008). Numerical models predict that, 
under some physical assumptions, merger remnants can con- 
tain extended gas disks (e.g. Hopkins et al. 2009). However, 
Bournaud & Elmegreen (2009) argued that the overall spatial 
distribution of the mass in high-redshift disks is inconsistent 
with the majority of them being merger remnants, suggesting 
that their mass assembly was dominated by relatively smooth 
mass infall along gas streams (see also Dekel et al. 2009b). 
This is because the morphology of these disk galaxies, in par- 
ticular their supermassive star-forming clumps, cannot be ac- 
counted for if the main assembly channel is by mergers. Our 
CO observations add further support to this picture. Indeed, 
we find extended CO reservoirs with spatial sizes consistent 
with what is measured from their stellar light. The end prod- 
uct of galaxy mergers are characterized by highly concen- 
trated molecular gas distributions, even when some molec- 
ular gas is found in the outer regions. This holds both for 
early-stage interactions (pre-merger galaxies like Arp 105 - 
Due et al. 1997), on-going mergers (Uke NGC 520, Yun & 
Hibbard 2001), and merger remnants (like NGC 7252, Hi- 
bbard et al. 1994). This dynamical effect results from the 
gravitational torquing of the gas, which is well reproduced by 
numerical simulations (Barnes & Hernquist 1991; Di Matteo 
et al. 2008) and is independent of the gas fraction. If ma- 
jor mergers played a role in the formation of z = 1.5 BzK 
galaxies, their large molecular disks must have been replen- 
ished by the smooth infall of baryons at a rate higher than 
the merger-driven mass assembly in order to keep extended 
gas disks. Cold flow-driven galaxy formation, indeed, does 
predict a large population of massive, extended star-forming 
disks at high redshift with high gas fractions (e.g., Dekel et al. 
2009b). 

7.4. Global properties of star formation at high redshift 

In a simple cartoon, the z = 1 .5 BzK galaxies appear to be 
disks like local spiral galaxies, but containing a substantially 
higher gas content and much larger star-forming clumps. A 
similar comparison holds for SMGs vs. local ULIRGs, with 
the former being the scaled-up version of the latter but hav- 
ing overall higher gas content (Tacconi et al. 2006). Most of 
the observed properties of the two classes can be understood 
within this scheme: there is a population of disk galaxies hav- 
ing low excitation CO, spatially extended gas reservoirs, low 
SEE (and hence long gas consumption timescales), low SSER 
and optically thin UV emission. These are the familiar spi- 
rals in the local Universe, and likely most of the near-IR (or 
optical) selected disk galaxies in the distant Universe. Con- 
versely, there are galaxy populations with opposite proper- 



ties: high CO excitation, compact gas reservoirs, high SFE, 
high SSER and with optically thick UV emission. These are 
the local ULIRGs and high redshift SMGs, both classes likely 
corresponding to major (wet) mergers of massive star forming 
galaxies. 

7.5. New constraints for galaxy formation models and the 

IMF 

Based on our results of a remarkably low scatter in observed 
SFE and gas fractions for the CO-detected BzK galaxies, we 
suggest that the molecular gas content plays an important role 
in driving the tight stellar mass-SER correlation in high red- 
shift sources. A simple scenario may be that the gas proper- 
ties are tightly connected to that of the underlying dark mat- 
ter halo, possibly via the regulation of cold gas accretion rates 
through cosmic times (e.g., Keres et al. 2005; see also Bouche 
et al. 2010). This in turn could regulate the amount of star for- 
mation in each halo and its time integral, i.e. the stellar mass, 
which, as a result, are tightly connected. Ultimately, there 
might be a stellar mass-SER correlation because the gas con- 
tent of galaxies tightly correlates with the masses of hosting 
dark matter halos over cosmological timescales. 

Galaxy formation models based on numerical hydrody- 
namic codes (e.g., Einlator et al. 2006) and semianalytical 
renditions (e.g., Kitzbichler & White 2007; see Fig. 17 and 18 
in Daddi et al. 2007) had indeed predicted the existence of 
tight correlations between stellar mass and SER before they 
were recovered by observations, with roughly a correct slope 
(or possibly slightly shallower, see, e.g., Daddi et al. 2007a). 
This is probably one of the rare cases where models of galaxy 
formation have preceded observations, showing truly predic- 
tive power Also, such models had qualitatively predicted an 
increase in the normalization of the correlation toward higher 
redshifts. 

However, existing models fail in reproducing the actual 
magnitude of the normalization increase to z ~ 2 (falling short 
by factors of 3-5 compared to observational estimates; see 
e.g. Daddi et al. 2007a; Dave 2008). The SERs estimated 
for distant galaxies are generally too high compared to the 
prediction of models, a fact that has been long known for 
SMGs (Bough et al. 2005). Dave (2008) suggest that this dis- 
crepancy between models and observations can be solved by 
adopting substantial evolution in the stellar initial mass func- 
tion (IMF). This conclusion was also reached by Baugh et 
al. (2005) in their modeling of SMGs (i.e., for galaxies with 
extreme SERs of typically 1000 Mp yr^^ that have space den- 
sities at least an order of magnitude lower than that of near-IR 
selected massive galaxies at z ^ 2). However, new theoreti- 
cal frameworks have recently been proposed in which feed- 
ing of star formation through cold flow filaments (Dekel et 
al. 2009a; see also Keres et al. 2005) can more easily jus- 
tify high star formation rates in distant galaxies. In addition, 
we have also found large gas reservoirs, reaching 10 "Mp in 
some cases, i.e. the gas accounts for a major fraction of the 
baryons in these sources. We are thus witnessing a phase in 
galaxy formation in which massive galaxies were truly gas- 
dominated systems. These high gas masses can explain the 
high observed SERs and one does not need to advocate non 
standard IMEs to reduce the SER, at least for the normal, near- 
IR selected (BzK) galaxies. We note that quantitative predic- 
tions of the models by Oppenheimer & Dave (2006; 2008) 
on the gas fractions in galaxies at z ^ 2 (with similar stellar 
masses to our CO-detected galaxies) result in expected gas 
mass fractions of order of 15%, i.e., much smaller than what 
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we derived. Therefore, our study suggests that the molecu- 
lar gas fraction in galaxy simulations should be increased in 
order to reproduce the observed IR luminosities, without nec- 
essarily introducing a top heavy IMF. 

8. SUMMARY 

We have presented results from a CO[2-l] survey of nor- 
mal, BzK-selected star-forming galaxies at z ~ 1.5. The anal- 
ysis of the integrated spectra and the spatially-resolved CO 
observations was made in combination with high-resolution 
optical imaging and multi-wavelength data. We also used 
appropriate numerical simulations of clumpy disk galaxies 
(BEE07) to derive dynamical masses and compare them with 
the observed CO luminosities. The main findings and impli- 
cations can be summarized as follows: 

• All observed galaxies were detected with high S/N ra- 
tios. We detected extended CO reservoirs with typical 
sizes of ^ 6-1 1 kpc, and presented evidence for rota- 
tion from spatially-resolved CO observations. This is 
indicative of the molecular gas being situated in disks 
that are similar in size and shape to the UV rest-frame 
morphology. 

• We have provided direct evidence for a high aco con- 
version factor (relating CO luminosity to gas mass), 
based on the dynamical masses of the systems. Our re- 
sults show that the conversion factor of BzK galaxies is 
high, similar to that of local spirals. Our derived aco is 
a factor of 4-5 larger than what has been derived for lo- 
cal ULlRGs and typically used for high redshift SMGs. 

• Our observations imply high gas fractions, ~ 50-65% 
of the baryons in our z ~ 1 .5 galaxy sample, and qui- 
escent star formation with relatively low SFE and long 
gas consumption timescales (~ 0.5 Gyr). 

• Our results also point to a strong correlation between 
the molecular gas content and other physical proper- 
ties of these sources. Once the stellar mass and SFR 
of a given galaxy are known, the CO luminosities can 
be predicted to within 40% rms. This and the over- 
all similar SFEs seen in these galaxies and local spirals 
leads us to propose a relation between CO luminosity 
and bolometric luminosity that could allow one to sys- 
tematically estimate the CO luminosity of massive disk 



galaxies, presumably at any redshift z ^ 2 as the rela- 
tion can fit local and z = 1 .5 observations with an over- 
all scatter of only 0.25 dex. 



• Our sample consists of normal star-forming galaxies at 
high redshift, and we have shown that it is represen- 
tative of the majority of z > 1 near-lR selected star- 
forming galaxies. Our results shed more light on the 
nature of these high-redshift disk galaxies. In particu- 
lar, the large spatial size of the molecular gas reservoirs 
suggest that they did not assemble mostly through past 
episodes of violent mergers. Instead, our observations 
are more consistent with smooth mass infall along cos- 
mic flows and subsequent internal evolution as the main 
drivers of disk galaxy formation and star formation at 
high-redshift. 
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